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In the holonomic approach to quantum computation information is encoded in a degenerate 
eigenspace of a parametric family of Hamiltonians and manipulated by the associated holonomic 
gates. These are realized in terms of the non-abelian Berry connection and are obtained by driving 
the control parameters along adiabatic loops. We show how it is possible, for a specific model, to 
explicitly determine the loops generating any desired logical gate, thus producing a universal set 
of unitary transformations. In a multi-partite system unitary transformations can be implemented 
efficiently by sequences of local holonomic gates. Moreover a conceptual scheme for obtaining the 
required Hamiltonian family, based on frequently repeated pulses, is discussed, together with a 
possible process whereby the initial state can be prepared and the final one can be measured. 
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The field of quantum information and computation 
(QC) Q brings together ideas and techniques from 
very different areas ranging from fundamental quantum 
physics to solid-state engineering and computer science. 
QC synergetically benefits from all these contributions 
and conversely quite often offers fresh viewpoints on old 
subjects. Recently it has been suggested Q that even 
tools related to gauge theories ||] might play a fruitful 
role in the arena of QC. Indeed in ref. || the possibility of 
realizing quantum information processing by using non- 
abelian Berry holonomies induced by moving along 
suitable loops in a control space M. has been analysed. 
The computational capability stems from the features of 
the global geometry of the bundle of eigenspaces associ- 
ated with a family T of Hamiltonians parametrized by 
points of Ai. The geometry is described by a non-trivial 
gauge potential A or connection, with values in the alge- 
bra u(n) of anti-hermitian matrices (n is the dimension 
of the computational space). Since the unitary transfor- 
mations realizing the computations are nothing but the 
holonomies associated with the connection A, this con- 
ceptual framework for QC is referred to as Holonomic 
Quantum Computation (HQC). In a sense HQC can 
be considered as the (continuous) differential-geometric 
counterpart of the (discrete) topological QC with anyons 
described in refs. 



In this paper we shall provide further analysis of this 
proposal. After concisely reviewing the conceptual basis 
of HQC, we shall show how, in a specific relevant model, 
one can explicitly determine the sequence of loops nec- 
essary for generating any given quantum gate. Then we 
shall introduce HQC models with a natural multi-partite 
structure and discuss how this bears on the question of 
complexity. Finally we shall discuss how in principle one 
can implement HQC by repeated pulse control of a sys- 
tem with degenerate spectrum. 



Let us begin by recalling the basic ideas of HQC [||. 
Quantum information is encoded in a n-fold degener- 
ate eigenspace C of a Hamiltonian Hq, with eigenvalue 
Eq. Operator, Hq, belongs to a family T = {H\}\ e _M, 
Hq = H\ , in which no energy level crossings occur as A 
ranges over M.. In the following we shall satisfy this lat- 
ter condition by assuming, for simplicity, that the Hamil- 
tonians H\ are isospectral (H\ = U(X) HoU(Xy). The 
A's represent the "control" parameters that one has to 
drive in order to manipulate the coding states £ C. In 
general the points of M. , from the physical point of view, 
can be thought of as describing external fields, such as 
electric or magnetic fields, or couplings between subsys- 
tems. Let C be a loop in the control manifold A4, with 
base point A , C: [0, 1] m> M,C(0) = C(l) = A . We 
assume that C is traveled along slowly with respect to 
the longest dynamical time scale involved: in this case 
the evolution is adiabatic i.e., no transitions among dif- 
ferent energy levels are induced. If \ip)i n £ C is the 
initial state, at the end of this control process one gets 
\i>)out — e l£oT T a{C) \tp)i n . The first factor here is just 
an overall dynamical phase and in the following it will 
be omitted; let us just mention that such a decoupling 
of the fast dynamical evolution opens new possibilities 
about coherent and error avoiding encoding ^ . The sec- 
ond contribution, the holonomy Ta(C) £ U(n), has a 
purely geometric origin and its appearance accounts for 
the non-triviality (curvature) of the bundle of eigenspaces 
over M.. By introducing the Wilczek-Zee connection M 



At'-= <V> a (A)|^-|^(A)>, 



(1) 



one finds Ta(C) = Pexp J c A Q, where P denotes path 
ordering. The set Hol(A) := {Ta(C)}c C U(n) is known 
as the holonomy group [||. In the case in which it coin- 
cides with the whole unitary group U (n) the connection 
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A is called irreducible. In it has been argued that for a 
large enough control manifold, the irreducible case is the 
generic one, therefore one can in principle implement any 
computation over the code C just resorting to this very 
special class of quantum evolutions. 

Quantum Gates. A workable HQC model which repre- 
sents a natural non-abelian generalization of the original 
Berry phase, for which explicit construction of the holo- 
nomic gates is possible, is now discussed. The model 
is worked out with some details in that it is extend- 
able to the more general case when A4 is a coset space. 
The features of the construction presented are twofold. 
On the one hand it fully exploits the loop composition 
structure at the basis of the holonomy group, showing 
a procedure whereby loops can be decomposed into 2- 
dimensional components, simple to deal with. On the 
other this topological construction overcomes the diffi- 
culties connected with the path ordering prescription. 

Let us consider the Hamiltonian Hq = eq\ii +l)(n+ 1| 
acting on the state-space H = C™ +1 = span{|a)}™ij. 
We shall take as the family T the whole orbit G(Hq) := 
{U H W I U G U(n+1)} of H under the (adjoint) action 
of the unitary group U(n + 1). This orbit is isomorphic 
to the n-dimensional complex projective space: 



O(H ) 



U(n + 1) 
U{n) x U(l) 



SU(n + l) 
U(n) 



CP' 



The points of CP™ can be parametrized by the uni- 
tary matrices U(z) = Ui(zi)Ui(zi)...U n (z n ) where 
U a (z a ) = exp[G a (z a )] with G a (z a ) — z a \a)(n + 1| — 
z a \n + l)(ct\, z a = 9 a e l ^ a . The eigenstates of the 
rotated Hamiltonians are \a(9, </>)} :— U(9,(j))\a) — 
cos9 a \a) - exp(-icf) a ) sm9 a YJjta exp(i<j>j) s'mdj 
U. j>7>a ^s0 7 \j) and \n + 1(9, 4>)) := U(9,cf>)\n + 1) = 

E™=i ex P(*^) sin ^ ri 7 <j cos6 '7b'} J where 9 n+1 := tt/2 
and 4> n +i '= 0. Notice that for n = 1 the standard 
2-level model with (abelian) Berry phase is recovered. 
By using equation (Q) the components of the connec- 
tion can be now explicitly computed. The only non-zero 
elements of the matrix A e > 3 ((3 — 1, . . . ,n) are AJ^ for 
a = 1,. .. ,(3 - 1, given by 



A p 

A a0 



fi>~t>a 



(2) 



as well as A- n — —A/-. The antihermitian matrix A^p 

ctp pet 

has non-zero elements for a ~ (3 and a > a given by 



A't'P — _Aj(<t>a-<t>li) 



sm tin sm 



sm a & 



n < 

/3>7>a 



with n,s>-Y>3 cos^-y = 1, and for (3 > a and a > a by 



A VI3 — if 
A aa — le 



/3>7>/3 T 



sin 9 a sin 9 & sin 2 9p J^J cos # 7 J^J cos # 7 . 

0>j>a /3>7>a 



The A 8 ? 's and A^f> 's are the 2 n components of the un- 
valued connection over CP n . 

For generating a given quantum gate g G U (n) one 
has to determine a loop C g in Ai such that r^(C ff ) = g. 
Due to the non-abelian character of the connection such 
an inverse problem is in general hard to solve. To tackle 
it we choose specific families of loops {Ci}, that gener- 
ate holonomies from which one can eventually construct 
any U(n) transformation. To this end we consider the 
2-dimcnsional submanifolds in the 2 n-dimensional space 
(0,4>), spanned by two variables, {0p,(j)p) or (9p,9^), for 
specific values of j3 and /3. For these loops the line integral 
is given by jL A = § c (A 6 t 3 d9p + A x ?d\p), where A = 
or (j). From (0) we see that we can always choose the pa- 
rameters which define the position of the plane (9/3, Xs), 
where loop C lies, in such a way that the matrix A 6 f 
is identically zero. If one takes 9i — 0, Wi ^ (3, (3, ma- 
trices A S/3 and A 9 commute, so that we can calculate 
the integral and exponentiate avoiding the path ordering 
problem. 

In this framework, it is possible to identify first four 
families of loops in such a way as to produce the basis of 
four matrices (the Pauli matrices and the identity) of all 
possible two-by-two submatrices belonging in the alge- 
bra of U(2). The first choice is (9/3, (pp), where the non- 
zero component of the connection is A = —ism 2 9 p. 
The second choice is the loop on submanifold (9/3, <Ag) 
for (3 > (3, with 0b = tt/2, giving a different connec- 
tion with two non-zero elements, At'i — i sin 2 9a and 



73/3 



A 



1 = 

4(8 



—i. Of course the latter element will give zero 



when integrated along a loop. For (3 < f3 both matri- 
ces are identically zero, and give rise to trivial holonomy. 
With these two connections and for appropriate loops 
one can obtain all possible U(n) diagonal transforma- 
tions. For loop Ci G (00, <f>p), T A (C X ) = expHljSXfllii], 
Ei denoting the area enclosed by C\, on the S 2 sphere 
with coordinates (20p,<j)p). For Ci G (0/3,(j>g), T A (Ci) = 
exp[i|/3)(/?|E 2 ]. Recalling the constraint (3 > (3, we see 
that one can produce n — 1 distinct holonomies from Ci . 

To obtain the non-diagonal transformations one has 
to consider a loop on the (00,0s) plane, with 0i = 
for all i ^ (3, [3. Then the only non-vanishing elements 
of the connection are A% = e^" - ^ sin0 fl = -A%. 

(3/3 p (3/3 

By choosing further the (9/3, 9/3) plane at = 0^ = 
the holonomy becomes, for C3 G (0/3,00)(i> ? =<j>s=o, 
T A (C 3 ) = exp[-(|/3}(/3| - |/3)(/3|)E 3 ], while at fa = tt/2 
and <j)p = 0, for C 4 G (0/3, 9p) 4>l3=7r/2 ^ j3=G , T A (C 4 ) = 
cxp[— i(\(3)((3\ + \f3)((3\)T,4], where E is the area on the 
sphere with coordinates (tt/2 — 00,0s). Note that any 
loop C on the (90,Xs) plane with the same enclosed 
area (when mapped on the appropriate sphere) Ec will 
give the same holonomy independent of its position and 
shape. These four holonomies restricted each time to a 
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specific 2x2 submatrix, generate all U(2) transforma- 
tions. Finally it is easy to check that in this way one can 
indeed obtain U = exp[//j Tj], where Tj (J = 1, . . . , n ) is 
a u(n) generator and \Xj an arbitrary real number. There- 
fore any element of U(n) can be obtained by controlling 
the 2n parameters labelling the points of CP™. 

It is instructive to consider the form that the Hamil- 
tonian family T , takes when restricted to the particu- 
lar 2-submanifolds. For the loop C\ (similarly for G-x) 
one finds H x = -e /2B(2^,^) • a for B(Q h 4>j) = 
(sin Oi cos (j>j, sniOi Bvafy, cos 9i) T , where the only non- 
zero elements are on the (3-th and (n + l)-th row and 
column. Hi generates an abelian CP 1 phase between the 
states |/3) and \n+ 1). On the other hand for the path C3 



(similarly for C4) we have H 3 = eoB(9i3,8p)B 



op) 1 ;, 

where the non-zero elements connect the states |/3), \(3) 
and |n + 1). In this Hamiltonian there is direct coupling 
between three states, giving rise to a non-abelian inter- 
action. 

In order to represent a 2-qubit system we have to con- 
sider the control manifold CP 4 . The holonomies in this 
case are 4x4 matrices, and we take as a representa- 
tion basis of the unitary transformations the qubit basis 
|00), 1 1 ) , 1 10) and |11). From the general scheme above 
it follows that by appropriate control of the parameters 
(8, (f>) for obtaining various loops C, we can generate all 
possible £7(4) rotations i.e., any logical gate, in partic- 
ular single-qubit rotations, and two-qubit gates such as 
the controlled operations XOR and CROT. (Uxor '■— 

|0)(0|®1 + |1)(1|®C7 X , Ucbot := |0)(0|®1+|1)(1|®<7,). 

For single qubit rotations we consider three unitaries 
U a = T A (C a ), (a = 1,2,3) with Ci G (0i, </>i)e 2 =0 2 =o, 
C 2 G (6*i, 02)e 2 =7r/2.0 1 =o, C3 G (#i,6» 2 )0 1= 2= o- For areas 
£ 2 = Si let U ph i = (U^U^U^)- 1 , equal to 



Uphi — 



COS £3 



sin£3e- 2iSl 



smS^e 2 ^ 1 cos £ 3 




(3) 



With loops on the (#3, <fe)e 4 =0 4 =o, (^3, <M0 4 =7r/2,0 3 =O) 
(#3, #4)03=04=0 planes and spanning the same areas as 
before, one can produce the rotation U p h2 on the lower- 
right block, and with their product obtain the phase ro- 
tation of one qubit 1 ® U q = U p h\U v h2- In order to 
perform the phase rotation of the other qubit all we 
need is the swapping operator, S, acting on the two 
qubits by S \ip) ® \4>) = \<j)) ® \ip), S is given in terms 



of holonomies, by S — r^C^)!^ 



=71-/2 



t/2' 



where C 4 G (6> 2 , #3)02=7.72,03=0 and C 2 G (£>i, <h)e 3 =w/2 
and <^3)e 3 =7r/2- Hence, the phase rotation of 
the other qubit is given by U q ® 1 = 5(1 ® 
Uq)S. The controlled rotation gate CROT, is given 
by Ucrot = r A (Ci)| Si=w/2 r A (C 2 )| S2=7r/2 , where 
Ci G (6> 4 ,<M and C 2 G (6*i, 4 )e 4=Tr / 2 . In ad- 
dition the 'exclusive or' gate XOR can be real- 
ized as Uxor = ^a{Ga)\±^/2 r A(C2)| S2=7r / 2 , where 



C 4 G (6*3, ^4)03=^/2, 04=0 and C 2 G (8i,<j> 3 )e 3 =w/2 and 
(6*1, 04)e 4 =Tr/2- 

Complexity. So far the coding subspaces analysed for 
HQC do not necessarily involve quantum entanglement; 
TL could be the state-space of a single quantum system. 
This is due to the fact that C does not have a built-in 
tensor product structure, thus, in general, it cannot be 
naturally interpreted as the state-space of a multi partite 
system. This latter feature, however, is one of the essen- 
tial ingredients that make QC more efficient than clas- 
sical computation. Indeed, from the above construction 
for the CP n -model it is easy to realize that the num- 
ber of elementary loops in 2-submanifolds necessary for 
implementing single and two-qubit operations scales ex- 
ponentially as a function of the number log 2 n of encoded 
qubits. To overcome such a difficulty one has simply to 
consider an Hamiltonian family, acting on the state-space 
of a multi-partite quantum system [not just isomorphic 
to], with a special structure allowing for local QC's to be 
holonomically performed. Then global QC's involving 
non-trivial actions over many qubits, can be efficiently 
realized in the standard way Q . One possible formaliza- 
tion of this idea is contained in the following theorem: 

Let H := <8>™=iWj ® H a , H 2 = C d , and H a = 
span{|— ), |+)} a single-qubit ancillary space. We set 
H a := ea z G End(ft Q ) @. Let moreover H(X) := 
'}2 i< jVij(\ij), where the Ay's belong to local control 
manifolds Mij. Suppose that Vij(X) G End(Hi ® Hj <8> 
H a ), and VijiO) = H a + H t + Hj, where Hj G End(W,-) 
is such that Hj \a)j = 0, (a — 0, 1) and that the family 
{Vij(X)}\eMij allows for universal HQC over the degen- 
erate eigenspaces Cfj :— span{|a)i ® \f3)j ® |±) : a,/? = 
0, 1} S C 2 ®C 2 . Then the family {H(X) /AG U t<j Mij} 
allows for efficient universal HQC over the n-qubit codes 
C± := span{®5L 1 |a i > ® |±> /a, = 0, 1} S (C 2 )®". 

The proof of the latter proposition proceeds as follows. 
First one observes that the subspaces C^j are degener- 
ate eigenspaces of Vy (0) corresponding to the eigenval- 
ues ±£. Hence the C s are 2™-dimensional eigenspaces 
of H(0) with eigenvalues ±n(n — 1) e/2. By assump- 
tion for any pair (i, j) of subsystems pT[ one can gen- 
erate any unitary transformation over C~j by adiabatic 
loops 7^ in A4ij. Keeping all the remaining A's at 0, one 
has a trivial action over the other factors of C , while 

T A hi 3 )-- ®k\a.k) ® |±) l-> Yspipj Up i /3 j ,a iaj ('yij)®k |0!fc)® 

|±), where ctk — ctk for k ^ i,j and ctk = Pk when 
k = In particular one can obtain a universal set of 
gates e.g., XOR^, and single-qubit operations, by using 
a fixed amount of resources. The claim then follows by 
well-known universality results for QC |l2]. The above 
scheme involves the use of an ancilla and requires con- 
trollability of three-body interactions, extremely difficult 
to achieve in practice. In this respect a simplification, 
involving just two-body interactions, can be obtained by 
considering N subsystems with d-levels O] . 
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Implementation. Now we discuss how one could in 
principle implement the holonomic loops even when the 
parametric Hamiltonian family T (or a part of it) is not 
available from the outset. We shall resort to ideas of 
quantum control theory in a way quite similar to the one 
adopted for symmetrization procedures JH] and decoher- 
ence control in open quantum systems JL5|. Suppose that 
an experimenter has at disposal the following resources: 
i) A quantum system characterized by the Hamiltonian 
Hq admitting a n-fold degenerate eigenspace C; ii) The 
way to turn on and off a set of interactions very quickly 
(with respect to the time-scales associated with Hq) in 
such a way that a set of unitary "kicks" K :— {U\}\ eM 
can be realized. Let T = N At and to — 0, ij+i = 
ti + A t (i = 1, . . . N — 1) be a partition of the time inter- 
val [0, T\. Now let the system evolution to be as follows: 
at any time ti the experimenter kicks the system with 
the pulse U- +1 Ui where Ui := C7(Aj) is a unitary chosen 
from the set K (Uq = Un = !)• Between the kicks the 
system evolution is unperturbed, U{At) = e~ lH ° At . The 
global evolution is then given by 



N-l 



U NAt (T)=T l[ UiU(At)U}, 



(4) 



where T denotes time-ordering. By considering the limit 
A T ^ 0, N i-> oo, with N At = T one gets U NAt (T) -> 
Texp -i J T dtH(t) where H{t) := U (X(t)) H U (X(t)Y . 
In particular by making the function X(t) vary adia- 
batically, one can obtain the desired holonomic evolu- 
tion. This scheme is based on a strong separation be- 
tween time-scales: each U\ has to be enacted impul- 
sively whereas the characteristic variation time of the 
control parameters A has to be slow enough to satisfy 
the adiabaticity requirement. More precisely, if Tk de- 
notes the kicking time, us the (highest) frequency associ- 
ated to the dynamics generated by Hq and t\ the time- 
scale over which the function \(t) varies, one must have 
Tfc < At <C us^ 1 <C t\. Notice that the pulses U\, are 
not required to be a universal set of gates for QC; here 
they represent an extra resource needed for implementing 
HQC when T is missing. 

Finally let us briefly consider the problem of codeword 
preparation and measurement. In order to encode the 
initial state into the degenerate subspace C, or to make 
the measurement on the final state, it would be useful to 
lift the degeneracy. Indeed in this way one would be able 
to distinguish energetically the different coding states. 
This characteristic is quite often desirable from the ex- 
perimental point of view in that one can resort to proce- 
dures involving energy transitions with state-dependent 
frequencies. The idea is to lift the degeneracy between 
the coding states, for example by switching on coherently 
an external (generic) perturbation. The basis states \ip a ) 
of C are mapped onto a set of states \ip' a ) that are no 
longer energy degenerate. Preparation/measurement are 



then performed and eventually degeneracy is coherently 
restored by switching off the perturbation. 

Summary. In this paper we have provided further anal- 
ysis of the proposal for Holonomic Quantum Computa- 
tion of ref . 1 1] . We have explicitly designed control- loops 
whose holonomies generate universal gates for a CP™ 
control parameter manifold. The basic idea is to as- 
sociate to the U (n)-generators computable transforma- 
tions obtained by loops on 2-dimensional subspaces of 
the control manifold. Explicit realization of two qubit 
gates has been given with the indication of which partic- 
ular loops the experimenter has to perform. In terms of 
such elementary holonomic gates we analysed the com- 
plexity problem and we showed how to achieve efficient 
implementation of quantum computing by resorting to 
a HQC model involving only local interactions. Some 
implementative issues have been addressed, and we de- 
vised a scheme based on repeated pulses for realizing the 
parametric family of isospectral Hamiltonians required 
for HQC. Finally, we briefly indicated how to prepare 
the initial state and how to measure the final one by co- 
herently switching on and off the energetic degeneracy of 
the computational subspace. 
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